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Abstract 

The motion of a Schwarzschild black hole with velocity vq = A)C through a constant magnetic 
field Bq in vacuum induces a component of the electric field along the magnetic field, generating a 
non-zero second Poincare electromagnetic invariant *F-F ^ 0. This will produce (e.g., via radiative 
effects and vacuum breakdown) an electric charge density of the order of p m a = Bo/3o/(2ireRG), 
where Rq = 2GM/c 2 is the Schwarzschild radius and M is the mass of the black hole; the charge 
density p- m & is similar to the Goldreich-Julian density. 

The magnetospheres of moving black holes resemble in many respects the magnetospheres of 
rotationally-powered pulsars, with pair formation fronts and outer gaps, where the sign of the 
induced charge changes. As a result, the black hole will generate bipolar electromagnetic jets each 
consisting of two counter-aligned current flows (four current flows total), each carrying an electric 
current of the order / « eBoRcPo- The electromagnetic power of the jets is L rj (GM) 2 B 2 j3 2 /c 3 ; 
for a particular case of merging black holes the resulting Poynting power is L ps (GM) 3 Bq/(c 5 R), 
where R is the radius of the orbit. 

In addition, in limited regions near the horizon the first electromagnetic invariant changes sign, 
so that the induced electric field becomes larger than the magnetic field, E > B. As a result, there 
will be local dissipation of the magnetic field close to the horizon, within a region with the radial 
extent AR ~ RgPo- 

The total energy loss from a system of merging black holes is a sum of two components with 
similar powers, one due to the rotation of space-time within the orbit, driven by the non-zero 
angular momentum in the system, and the other due to the linear motion of the black holes 
through the magnetic field. 

Since the resulting electrodynamics is in many respects similar to pulsars, merging black holes 
may generate coherent radio and high energy emission beamed approximately along the orbital 
normal. In addition, merging black holes may produce observable wind-driven cavities. 
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PACS numbers: 04.30.Tv, 95.85. Sz 



I. INTRODUCTION 

Observations of an electromagnetic signal accompanying black hole merger is most de- 
sirable, as it will provide crucial information on the location and the physical properties of 
the event. Two types of electromagnetic signal can be expected from merging black holes. 
First, merging black holes induce perturbations in the surrounding gas [1-4]. The resulting 
electromagnetic signal is then subject to great uncertainty and naturally depends on the 
complicated non-linear fluid behavior of the system. One of the problems is that at late 
stages of the merger there should be little gas inside the orbit since the timescale for shrink- 
age of the binary orbit by gravitational wave radiation becomes shorter than the timescale 
for mass inflow due to viscose stresses in the disk [3]. It is then hard to excite transient 
dissipative processes in a faraway accretion disk. 

Alternatively, external gas can support electric currents that create large scales magnetic 
fields. Motion of black holes in this externally supplied magnetic field can then lead to 
electromagnetic extraction of energy. Qualitatively, there are two distinct possibilities for 
the electromagnetic extraction of energy from spiraling black holes. First, the system of two 
orbiting black holes possesses non-zero angular momentum, which induces rotation of space 
time. Rotating space-time can generate electromagnetic outflows, in a manner similar to 
the classical Faraday disk. This is the physics behind the Blandford & Znajek [5] process 
of extracting the rotational power of a black hole. Below we refer to this mechanism as the 
Faraday disk mechanism. This mechanism has been previously considered in Ref. [6]. For 
a given magnetic field, the resulting power can be estimated using the Faraday disk scaling, 
Lem,f ~ R 2 B 2 c(QR/c) 2 , where R is the orbital radius and Q is the typical angular velocity 
of the rotation of the space-time within the black holes' orbit, Q ~ (GM) 3 ^ 2 /(c 2 R 5 ^ 2 ). The 
resulting Poynting flux is then 

g 3 m 3 d2 

Lem,f ~ B (1) 

In a separate, physically distinct process, which we consider in this paper, the linear 
motion of a Schwarzschild black hole will generate the electric potential drop across the 
black hole, so that a black hole will effectively operate as a unipolar inductor, in a way 
somewhat similar to the planet Io moving in Jupiter's magnetic field [7]. Classically, the 
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motion of a conductor though the magnetic field generates in the frame of a conductor an 
electric field E = — v x B. This induced electric field will generally have a normal component 
to the surface of the conductor. As a result, surface charges will be generated; they will 
produce their own electric field, now with a component parallel to the initial magnetic field. 
This electric field drives currents along magnetic field lines; dissipation of these currents is 
responsible for non-thermal radio through X-ray emission of the Jupiter magnetosphere. An 
important difference of the unipolar induction mechanism considered in the present paper 
from the classical unipolar inductor is that in the case of a black hole no physical charges 
are needed to produce E\\. Parallel electric field appears in complete vacuum due to the 
curvature of space. 

A qualitative estimate of the resulting Poynting power may be obtained from the following 
reasoning. For a conductor of length / moving with velocity Vq = /3qC through magnetic field 
B the resulting potential difference A$ ~ (3 LB. If the resulting outflow is relativistic, the 
electromagnetic power can be estimated as L E m ~ A$ 2 c = /3qBq1 2 c (see also [8]). In the 



case of orbiting black hole, estimating / ~ R G (the Schwarzschild radius) and (3 = yRc/R, 
the expected power of unipolar inductor turns out to be the same as that of the Faraday 
disk mechanism, Eq. (1). 

Thus, the estimates of the electromagnetic powers due to rotation of the space-time and 
due to linear motion of a black hole with Keplerian velocity turn out to be similar, given by 
Eq. (1); we view this as a coincidence. Though in both cases the power eventually comes 
from the inductive electric field, the underlying physics is different in the two case. One 
mechanism requires non-zero angular momentum, while the other does not. A total energy 
loss from a system of merging black holes is a sum of two components with similar powers, 
one due to the rotation of space-time within the orbit, another due to linear motion of the 
black holes through magnetic field. 

Previously, in Refs. [9-11] a number of force-free simulations of black hole magnetospheres 
were performed. In the case of the orbiting black holes the authors mostly studied the 
electromagnetic power due to linear motion of black hole, and not due to rotation of the 
space-time within the orbit. The present paper offers explanations and parameter scalings 
of these numerical simulation. 
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II. STATIC ELECTROMAGNETIC FIELDS IN SCHWARZSCHILD METRIC 



Consider vacuum stationary homogeneous electromagnetic fields in Schwarzschild metric. 
Though the structure of a constant electromagnetic fields in Schwarzschild metric is well 
known [12-15], here we briefly re-derive it here for completeness. Adopting a Schwarzschild 
metric 

ds 2 = -a 2 dt 2 + \dr 2 + r 2 ( d6 2 + sin 2 Odcf) 2 ) (2) 
where a = \Jl — 2M/r, the relevant vacuum Maxwell equations 

% (V=gF^) = o, (3) 

(F^u = A[^ u } is Maxwell tensor) for the non-vanishing components of the vector potential 
A (r,9') and A^r, 9) (here 9' and 9 are angles with respect to the axes aligned with the 
electric and magnetic field at infinity) give 

a 2 d r (r 2 d r A ) + -J—de,(sm9'd e ,A ) = 
sin 9' 

r 2 d r a 2 d r A^ + sin 9d e (^—d e A^ =0. (4) 

(In vacuum, the same relations hold for the dual Maxwell tensor *F flu ; in that case the 
equations for A and A^ switch.) The potentials corresponding to constant fields at infinity 
are 

A = E {r-2M) cos 9' 

A<t, = ^r 2 sin 2 9 (5) 

The electromagnetic fields 

E = {I /a) V A) = £ (cos 9'e r - a sin 9'e e ) 

B = V x A^e^ = B (cos 9e r — a sin 9e e ) (6) 

where V is a covariant derivative, with a corresponding unit radial vector e r = ad r . The 
structure of the electric and magnetic fields is the same, as follows from the duality trans- 
formation in vacuum. 

The above relations can be obtained in a more conventional way by using the alternative 
3 + 1 formulation of the General Relativity [16], in which case the Maxwell equations in the 



more general Kerr metric take the form 

V • E = 4np 
V-B = 

V x (aB) = 4?raj + D t E 

V x (aE) = -AB (7) 

where D t = d t — £g is the total time derivative, including Lie derivative along the velocity 
of the zero angular momentum observers (ZAMOs). For Schwarzschild metric /3 — and 
in the stationary case D t = 0. Electromagnetic fields in Eq. (7) are those measured by a 
stationary local observer in terms of local time. The fields aE and aB are those measured 
by a local observer in terms of Schwarzschild time t. 

Since we expect that the non-zero charge will eventually be created, we give here the 
complete Laplace equation for the electric potential A in presence of non-zero charge density 

V-VA = 4vrp 
a 

^dri^drAo) + - 1 d {sm OdgAo) + - 2 \ 2 dlAo = 4vr^ (8) 
r z a z r sin 9' a l r l sin 9 v a 

cf. [17, Eq. 10]. 

The shape of the field lines (6) is given by 

sin6 = 2e-^ — r -^- (9) 

where r_i_ i0 is the radial cylindrical distance from the axis of a given field line at infinity. 
The last field line that intersects the black hole at 9 = n/2 has initial r_|_ i0 = e 2 M/2 and 
is given by sin# = Note, that though the surfaces of constant magnetic flux are 

cylinders rsin^ = constant (corresponding to = constant), on the horizon the magnetic 
field in Schwarzschild coordinates becomes radial. This is naturally impossible at the point 
r = 2M,9 = 7r/2, but at that point the magnetic field is zero. This seeming inconsistency 
can be resolved in terms of embedding diagrams [13]. 

By construction fields (6) correspond to V • E = 0, so any local observer would measures 
zero charge density. On the other hand, judging by the shape of electric field lines in 
coordinates r — 9, the observer at infinity will infer a dipolar-like charge distribution 

cos 9' 

p eff = V nc -E = 2(l-a)— E (10) 



Here V nc denotes the non-covariant differentiation with respect to coordinates r — 9. The 
effective charge is concentrated near the black hole so at large distances the black hole 
effectively has a surface charge 

1 f°° \n(e/2) 
<j eff = — / Peff dr = -^^E cos 9' (11) 

47T J 2M 71 

We stress that there are no physical charges present: the shape of electric field lines in the 
chosen metric is modified by gravity, not electric charges. Also, the shape of field lines and 
the value of the effective surface density depends on the choice of coordinates and a given 
set of observers. 



III. BLACK HOLE IN ELECTROMAGNETIC FIELDS ORTHOGONAL AT IN- 
FINITY. 

Let us next assume that in a particular reference frame at large distances from the black 
hole the magnetic field is along — e z direction, while in this reference frame the electric field 
is zero. A black hole is moving orthogonally to magnetic field with velocity (3 = v /c along 
y direction. In the frame of the black hole, at large distances r>M, there is then a static 
electric field E = — v x B = /3 B e x , E = /3oB . (In this section we define B as the value 
of magnetic field in the frame where black hole is at rest; it is related to the value in the 
frame where electric field is vanishing by a simple Lorentz transformation). Thus, in the 
frame of the hole the electric and magnetic fields are given by Eq. (6), where 9' is a polar 
angle with respect to the x axis, cos 9' = sin 9 cos 0. The electromagnetic invariants are 

B 2 -E 2 = Bl(l - (3l) + 2 (/3 2 (cos 2 9 cos 2 + sin 2 0) - sin 2 O) ^ 

E B = - Det J^H = -cos0sin2^ OJ B 2 — (12) 

while the parallel electric field is 

M 

E\\ = - cos <f) sin 29 — - =£p (13) 

ry/l -2 sin 2 9M/r 

Eqns. (12) highlight two important points. The second Poincare electromagnetic in- 
variants is generally non-zero, E • B ^ 0. The first Poincare invariant changes sign on the 

2M (sin 2 - ft 2 (cos 2 cos 2 <p - sin 2 <£)) 

r = T=M (14) 




FIG. 1. Shape of magnetic (solid lines aligning with z axis) and electric (dashed lines, aligning 
with x axis) field lines in the x — z plane for vacuum Schwarzschild black hole. Red curve (in the 
online color image) starting at the point 6 = vr/2, r = 2M is the shape of the outer gap, where the 
sign of the induced charges changes. 

This occurs close to the horizon around points {6 ~ vr/2,0 f» 0,7r} within a region ^ w 
ffi/(l-0%), see Fig. 2. 

We expect that in astrophysical environment En will be decreased by pair creation (see 
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FIG. 2. 3D view of the black hole magnetosphere. Magnetic field (thick solid lines together with 
magnetic surfaces) is along z axis at large distance, the black hole is moving along the y direction. 
The central sphere of radius r = 2M is the Schwarzschild radius. Two limited regions near the 
equatorial plane approximately aligned with the x axis are bounded by the surfaces E = B. On 
the nearly conical surfaces originating at the magnetic equator the induced charge density is zero. 
This particular plot is for (3q = 1/2. 

§IV). Still, the perpendicular component of the electric field E± becomes larger than the 
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magnetic field at points 



2M sm 2 9\ 2 



fta 1 - 



2M sin 2 9 sin 2 



(15) 



The region bounded by this surface has a shape similar to the one where E = B, Eq. (14), 
see Fig. 3. 



y/M 




FIG. 3. 3D view of the surfaces bounding the regions where E > B {Left Panel) and E± > B 
{Right Panel) for /3q = 0.5. The central sphere of radius r = 2M is the Schwarzschild radius. False 
colors (in the online version) in this figure and Figs. 4 and 9 are chosen for presentation purpose 
only. 



Since we expect that a plasma will be generated due to vacuum breakdown, we can define 
the plasma drift velocity (assuming (3d < I, see below) 

(3d = E x B IB 2 = 9 - {a sin 9 sin (pe r + cos 6 sin (f>e e + cos ^e^) (3 (16) 

1 — 2Msin 9/r 

It is directed along y axis at large distances, but becomes complicated near the hole. Most 
importantly, the value of v is not bound. Close to the points 9 = ir/2, <p = 0, tt, the velocity 
diverges as (3d ~ Po/at- This is related to the previously mentioned fact that the second 
electromagnetic invariant changes sign close to these points. 
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IV. MAGNETOSPHERES OF MOVING SCHWARZSCHILD BLACK HOLES 



A. Induced charge density 

In the previous section we showed that the motion of a Schwarzschild black hole through 
magnetic field in vacuum generates non-zero second Poincare electromagnetic invariant (par- 
allel electric field) and regions where the first Poincare changes sign, E > B. In vacuum this 
would have been the end of the story, so to say. In reality, astrophysical plasmas tolerate 
neither large parallel electric field, nor regions of E > B: any stray particle will be acceler- 
ated by electric field and will produce an electron-positron pair via various radiative effects. 
The resulting pair plasma will try to screen the initial parallel electric field, by producing a 
charge density required to have Ey = 0. 

Thus, we expect that around astrophysical black holes moving in external magnetic field 
a non-zero charge density will appear, that satisfies the equation (8) with induced electric 
field 

E ind = £,|b (17) 



where 

1 



[— cos9e r + asinOeg) (18) 



yl - (2M/r) sin 2 9 

is a unit vector along magnetic field. The corresponding charge density, see Fig. 4 



1 a (1 + (M/r- 3/2) sin 2 9) _ ME 



Pind = V • E = — v . - ' - - . ' , \ — J cos(/>sinfl — (19) 

1 4tt tt (1 - (2M/r)sin 2 #) 2 Y r 2 y ' 

The typical value of charge density p is 

Po ~ 4^M " 2nc (20) 

This reminds the Goldreich- Julian charge density [18] if we substitute v /R G — > £l e ff- 
At large distance from the black hole the charge density is 

(1 + 3 cos 29) cos sin 9M 2 . , 

Pind ~ ~ 2 Po (21) 

On the magnetic equator 9 = it/ 2 the charge density is 

Pind,6=ir/2 = 2 ^ COS (22) 



a r l 
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FIG. 4. Surfaces of constant charge density pi n d = +| Const | (Left Panel) and pi n d = — | Const | 
(Right Panel) for Const = 2 x lCT 2 po and /3 = 0.5. 

It diverges on the horizon. The charge density close to the horizon pi n d,h is 



It diverges near the equator. 

The point r = 2M, 9 = it/ 2 is a special one, as can be seen from the fact that the limit 
of Eq. (23) for 9 — > n/2 does not coincide with the limit r — > 2M of Eq. (22). This can be 
traced to the fact that on the one hand, the lines of force must cross horizon orthogonally 
[e.g., 13], yet on the other hand they must lay on cylinders r sin 9 = constant (see Eq. (6)). 
These two conditions are inconsistent at the point 9 = n/2, r — > 2M. Recall also that the 
magnetic field is zero at this point, Eq. (6), while experiencing a kink, clearly seen in the 
embedding diagram [13]. 

B. Plasma response 

Similarly to the case of rotating neutron stars, two alternative possibilities exist with 
respect to the dynamical response of the system to the required charge density. First, the 
overall plasma distribution can provide a static distribution of the required charge density 



cos 4> sin 8 



(23) 



Pind 



cos 2 9 
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and the system will overall be quiet, without strong electromagnetic outflows. In case of 
pulsars, this approach is advocated in Ref. [19]. This scenario is supported by the fact that 
direct PIC simulations of the pulsar magnetospheres so far failed to produced an outflow: 
the system indeed just relaxes to a nearly steady-state configuration [20]. It is expected 
that secular instabilities of the resulting charged configurations [e.g., diocotron instability, 
21] may eventually lead to the formation of the jets (Spitkovsky, priv. comm.). 

Alternatively (and this viewpoint is supported by majority of pulsar theorists), static 
charge configuration cannot be established on magnetic field lines connecting to infinity, 
resulting in the formation of the wind [18]. Fluid simulation of pulsar magnetospheres 
support this paradigm [22]. Since in our case all magnetic field lines are connected to 
infinity, large parallel electric fields will lead to plasma outflow, which would qualitatively 
resemble plasma outflow on the open lines of pulsar magnetospheres. 

We accept the paradigm that in the case of black holes, similar to the pulsar magneto- 
spheres, parallel electric fields will lead to vacuum breakdown, generation of primary beam 
and a dense secondary plasma. In this Section we investigate the resulting electromagnetic 
structure of the black hole magnetospheres. 

C. Pair formation front 

A black hole moving with velocity v = (3 c through magnetic field B creates a potential 
drop of the order A$ m PqTgBq. For example, for a black hole of mass M = 10 6 M & rriQ 
moving on the orbit £ ^> 1 times larger than the Schwarzschild radius, i? rb = C r G, m 
a given magnetic field, the resulting potential is $ ~ 5 x 10 13 V (j^J m^ 1 / 2 . This will 
result in a particle energy, which is typically much larger than the one required to break 
the vacuum. For example, in pulsar magnetospheres, the pulsar death line corresponds 
approximately to $ 10 9 V [e.g., 23, 24]. Though radiative processes and the photon 
fields around pulsars and around black holes can be substantially different, in case of black 
holes the available potential is many order of magnitude larger than the one corresponding 
to the pulsar death line. This should ensure the vacuum breakdown and the formation of 
the secondary pair plasma in black hole magnetospheres. 

For a given parallel electric field, Eq. (12) and a given shape of field lines, Eq. (9), we 
can calculated the total potential drop along a field line A$ = / E\\ds, where ds is taken 
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along a given field line. For magnetic field lines that intersect the black hole horizon (those 
that have r± :0 < e 2 M/2) the lower limit of integration is on the horizon, while for those field 
lines than miss the hole the lower limit of integration is at magnetic equator. After a particle 
gains sufficient energy eA$ cr ;t to break the vacuum via radiative effects, a dense flow of the 
secondary particles will screen the resulting parallel electric field. We can calculate the shape 
of this so called pair formation front (equipotential surface) by requiring J E\\ds = eA<3> crit , 
see Fig. (5). 

Near the central field line the pair formation front is located at large distances from the 
black hole since close the axis the electric field is small, E\\ ps 2E cos^^^, where r± is the 
distance from the axis. To find the shape of a pair formation at large cylindrical distances, 
we note that magnetic field are nearly straight, while the parallel electric field bcomes 
E\\ ps sin 2$ cos (f>(ME /r). A fixed potential drop is then achieved on surfaces satisfying 

eA$ t 

sin 9 PPF = a 7~ cnt - . (24) 
2M£ o cos0 

Formally, the pari formation front extends to infinity. 



D. Large scale jets: quadruple current flow 

Accepting the paradigm that the plasma will be streaming along magnetic field lines 
nearly with the speed of light, the charge density (19) will generate a current flow in the 
black hole magnetosphere. At large distances the sign of parallel electric field and of the 
charge density depends on the quantity sin 9 cos 0, Eqns. (12,19). Thus, two counter-aligned 
currents propagating along z axis will be generated (four total separate current), separated 
the y — z plane, Fig. 6. 

Each quadrant of space carries a total current 

/ « r 2 d(t) f 72 sin OdOp = * M 2 Po = GMP ° B ° (25) 
j-tt/2 jo 4 16 

This current will produce toroidal magnetic field (for cylindrical radii M <C r± z) 

Mr 2 

At a given cross-section z = constant, the toroidal magnetic field will roughly correspond 
to two equal counter-aligned current flows. 
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12 3 4 

x/M 

FIG. 5. Structure of the equipotential surface in the x — z plane (solid line). It is calculated by 
requiring that the potential drop along a given field lines equals some fiducial potential J Euds = 
eA<I> cr it. Dashed lines are magnetic field lines. 

The regions of space with two counter-aligned currents propagating along ±z direction 
(corresponding to semi-spaces x > and x < 0) will interact with each other. Since the 
currents are due to charge-separated flow, the interaction is both electric and magnetic. 
Oppositely charged currents have charge density per unit length A ~ M 2 p . The total 
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FIG. 6. Value and direction of parallel electric field En in the plane y = 0. Shadowing intensity 
is proportional to the absolute value of the electric field. The electric field has opposite direction 
in the four quadrants; as a result two counter-aligned current flows in each ±z direction will be 
generated. 



electric force per unit length will be F e pa X 2 /r = M 4 pl/r. This electrostatic force will be 
nearly balanced (but not completely) by the magneto-static repulsive force of two counter- 
aligned currents /, Eq. (25). (The repulsive Lorentz force in this case will be smaller than 
the attractive Coulomb force by 1 — v^/c <C 1, where t>& is the velocity of the particles from 
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the primary beam.) 



E. Outer gaps 



On some field lines the charge density changes along a field line, Fig. 7. By analogy with 
pulsar magnetospheres we will call the surfaces p m d = the "outer gaps". The condition 
Pind = is satisfied on surfaces given by 



2Msin 2 
3 sin 2 6-2 



(27) 



The outer gaps touch the horizon at the equator and extend within polar angles arcsin y 2/3 < 
6 < 7r — arcsin \j2/3. At large distances the outer gap becomes a conical surface with 
cos 26 = ±1/3, see Fig. 2. In addition, some magnetic field lines have the sign of the 



O 




-0.1 



-0.2 



FIG. 7. Charge density along magnetic field field lines starting at 9 = tt/2 at different az- 
imuthal angles cj) = — (3/4)ir, —ir, 0, tt/2 (top to bottom) at large r). The zero is at radii satisfying 

- — + 1 = 0, r = 2.260M. 



4y/l — 2M 

4M 2 e V? (^-|) 



( v / f=23f+ v ^) 
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induced charge density changing two times along the field lines (this occurs along a narrow 
annular bunch of magnetic field lines that cross the horizon close to the magnetic equator). 



V. ELECTROMAGNETIC SIGNALS 

A. Overall power 

What electromagnetic signal can be expected from Schwarzschild black holes acting as a 
unipolar inductors? In a magnetically dominated medium with the effective impedance of 
free space, Att/c, the currents (25) will result in energy loss 

Lem,u ~ M 2 El = M 2 Blft (28) 

The power (28) approximately equals B 2 r 2 c with toroidal magnetic field given by the esti- 
mate (26). Equivalently, it arises due to initial electric field and the toroidal magnetic field 
(26) induced by the currents. 

In particular, for a binary BH before merger, we can estimate the electric field as 

E ~ B ^jM/R (29) 

where R is the radius of the orbit. The resulting power is then 

M*Bl _ [GMfB 2 
L EM ,u - = c5R (30) 

where we reinstated the Newton's constant and the speed of light. This is an estimate of 
the power lost by a black hole moving through constant magnetic field via unipolar inductor 
mechanism. It turns out to be of the same order as L EM;F , Eq. (1), the power dissipated 
by the rotation of space with the black holes orbit via the Faraday disk mechanism [6]. 
Still, we stress that these powers comes from related but physically separate mechanisms: 
a Faradey-type electromagnetic outflow considered in [6] and the linearly moving unipolar 
inductor discussed here (see also [8]). The power (30) is typically much smaller than the 
power carried away by gravitational waves at the later stages of the merger [6]. 

B. Estimates of the magnetic field 

The general expression for the Poynting power (30) depends on the strength of the mag- 
netic field, which in case of merging astrophysical black holes must be supported by the 
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accretion disk and thus depend both on the microphysics of the disk dynamo [25] and the 
dynamical evolution of the binary black hole-accretion disk system [2, 3]. 

Since the Poynting power (30) is the same as in the case of Faradey-type electromagnetic 
outflow driven by the rotation of space-time, we can use various estimates of magnetic field 
and the merger dynamics discussed in Ref. [6]. Here we briefly restate the points. 

Two stages of disk dynamics may be identified. At early times the loss of energy via 
gravitational radiation is slow enough, so that due to viscous diffusion the inner edge of 
the disk will be located close to the orbital radius, -R d ps R. At later times, for R d < 
$, d GM/c 2 , £ d 40 — 100 [3, 26], the binary will decouple from the disk, undergoing a merger, 
while the inner edge of the disk remains fixed at R d - 

There are several estimates of magnetic field, which generally give similar values [6]. For 
example, magnetic field can be estimated assuming that a fraction t]e of the Eddington lu- 
minosity is carried by magnetic field, B 2 ps rj E L Edd / (cR d ), (where R d ~ R before decoupling 
and R ~ i d Ra after decoupling) 

B = ,f^fl „ 3 x lO'G^^ (£)-' (31) 
Thus, after decoupling magnetic field remains nearly constant within the orbit, 

Bt-SOOGm^yj^^l (32) 
The total power (30) in this case then becomes 

^ V Idd2 <VE& 2 ( 33 ) 



L EM {GMf 
^Edd c 6 RR d 



The peak electromagnetic power is typically a fairly small fraction of the Eddington lumi- 
nosity Lem ~ 10~ 5 — 10~ 3 LEdd and can realistically be observed at cosmological distances 
only for merger of very massive black holes, M > 1O 8 M . 

Finally, comparing the power of the Schwarzschild black hole as unipolar inductor, Eq. 
(1) with the Blandford-Znajek power in a given magnetic field, L B z ~ a 2 B 2 M 2 (a is the 
black hole spin parameter), gives 

L EM _ M (M) 



L BZ a 2 R 

Thus, for fast rotating black holes, a ~ 1, the power of the unipolar inductor is subdominant 
for R 3> M. (In addition, since during the merger the inner edge of the disk is at fairly 
large radii, > 40 — 100-Rg, it is expected that the magnetic field on the black hole is smaller 
in the case of a merger). 
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C. Expected emission 



Eqns (30-33) provide estimates of the total electromagnetic power produced by a black 
hole in a form of a Poynting flux due to unipolar induction mechanism. A fraction of this 
power will be dissipated and converted into the observed radiation. Next, we discuss possible 
emission signatures. By analogy with pulsars, we expect several types of electromagnetic sig- 
nal from black holes as unipolar inductors, approximately corresponding to magnetospheric 
and plerionic emission. In addition, in case of black holes we expect emission from regions 
with E > B. 

1. Magneto spheric-type emission 

As we discussed above, the magnetospheres of black holes moving in magnetic field re- 
semble in many respects the pulsar magnetospheres. Thus we may expect somewhat similar 
radiative signatures, though the details of the emission mechanisms proposed below must 
naturally be calculated independently for black hole magnetospheres. Rotationally-powered 
pulsars produce, generally speaking, two types of radiation: coherent radio emission and 
high energy X-ray through 7-ray emission. 

As we discussed in §IV E, the sign of the induced charge density changes as a function of 
distance from the hole along a set of file lines, Figs. 2,7. One might expect that this will 
lead to effects qualitatively similar to the ones occurring in the so-called outer gaps in pulsar 
magnetospheres, where the sign of the Goldreich- Julian density changes [27]. In particular, 
recent Fermi observations of pulsars show that high energy 7-ray emission is generated at 
the outer gaps [28]. Similarly, we can expect that merging black holes will produce high 
energy emission, which will be preferentially beamed along the normal to the orbital plane. 
One might expect that, similar to pulsars, the emitted high energy power may reach tens of 
percent of the total electromagnetic power, Eq. (30). 

The radiation physics of pulsar outer gaps is not well understood at the moment [29]. A 
particularly important difference between the pulsar and black hole magnetospheres consid- 
ered here is that in case of pulsars the curvature radiation (in addition to inverse Compton 
processes) is an important ingredient. Magnetic fields in the magnetospheres of black holes 
will generally have drastically different - larger - radii of curvature, that would make curva- 
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ture radiation unimportant (this cane be easily demonstrated following the estimates of the 
inverse curvature emission given below). 

Another, in addition to curvature emission, important ingredient in the dynamics of 
the pulsar outer gaps is the inverse Compton (IC) emission. It will provide the dominant 
radiative friction effect to the particles accelerated in the black hole magnetospheres. Let us 
next estimate its properties. First, balancing the radiation friction force ~ 7 2 e 4 ?7 p h/(mgC 3 ), 
where 7 is the Lorentz factor of a particle and C/ p h is the radiation field energy density, with 
the acceleration rate by the induced electric field, eE c, the terminal Lorentz factor becomes 



7=meC V^ (35) 

Estimating electric field as E = (3 B = ^ Rq / R orh B and relating the radiation field energy 
density to Eddington luminosity, = L B dd/ (^R^c), we find 

= / GMB,nW \ ^ = e x l0 , Bl n mre i< (36) 
V e ' A m v J 

where, we remind, R or h = C r G- 

We can then estimate the IC power produce by the leptons with typical density given by 
Eq. (20) within a typical volume w Rq, 

L IC « ^B 2R l c = 1040 ergs' 1 BlmlC 1 (37) 

The typical frequency of IC emission is given by up-scattering of thermal photons of energy 
emitted at the inner edge of the accretions disk with energy 

, _ / GMM \ 1/4 

k B T « — 38 

where cjsb is Stefan-Boltzmann constant, k B is Boltzmann constant and M is accretion rate. 
Scaling accretion rate M with Eddington luminosity, rj M Mc 2 = L Edd , where r\ M ~ 0.1 is the 
accretion efficiency, we find 

/ 5 \ 1 / 4 

^fe 3 ) ^°eWW (39) 
The corresponding IC photon energy is 

e « 7 T = - — - =— — — = = 800 GeVfi^e £ 7 »7-i (40) 
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Luminosities (37) at energies (40) can be detected in reasonable time (e.g., 100 hours of 
observations) by atmospheric Cherenkov telescopes like HESS and VERITAS only within 
the few Mpc. 

Similarly to the case of pulsars, merging BHs may also produce high brightness coherent 
radio emission. The parallel electric field (17) will produce a primary beam of density (20) 
moving with highly relativistic velocity. Pair production by the particles from primary beam 
will generate secondary plasma with much higher densities. Such momentum distribution of 
particles may result in various types of plasma instabilities and generation of coherent high 
brightness radio emission [e.g., 30, 31]. 



2. Plerionic-type emission 

Most of the power (30) will leave the black hole region in a form of relativistic highly 
magnetized wind. Even though the instantaneous power (30) is typically much smaller than 
the Eddington power corresponding to masses M, Eq. (33), the total released energy can 
be substantial as we demonstrate below. Most of the energy is released before decoupling. 
Indeed, after the decoupling magnetic field is constant within the orbit, Eq. (32), so that 
the power (30) integrated over the merging orbit with R (i? 4 - — G3 Jf 3f ) ^ after decoupling 
gives 

E d = VE J {G ^^dt « VE L Edd ^ « BlRl (41) 

Thus, after decoupling the black holes spiraling in a constant magnetic field dissipate ap- 
proximately the magnetic field energy within the volume of the decoupling radius. 
On the other hand the total energy released, before decoupling is 

„ -^Edd-Ro GMm p R 

E t ot ~ Ve = Ve — (42) 

where Rq is the initial radius where the model becomes applicable. Thus, most of the energy 
is released at large separations of the black holes. 

As an estimate of the initial radius Rq where the model is applicable, one can compare 
the pressure created by the outflow, L em / (R%c) with the thermal energy density in the 
surrounding medium, ps n cx m p cj, (n cx is number density and c s is the speed of sound in the 
surrounding medium). This gives 

/ (CM) 4 \ 1/5 3 4/5 _ 2 /5 -1/5 1/5 
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JP ({GM) 9 ml\ 47 9/5 -2/5 6/5 -1/5 ,, n 

E t°t~ ^Tl 6 = 10 er S m 6 C S ,6 VE,-l n ( 43 ) 



1/5 



This is a substantial amount of energy, especially for mergers of more massive black holes. 
It can in principal be observed both via direct emission and by appearance of dynamical 
morphological feature resembling a wind-blown nebular in the central parts of galaxies. 



3. Emission from from regions with E > B 

Violation of the condition B > E implies large electric field in plasma that cannot be 
reduced by the relative motion of charges (recall that the drift velocity is independent of 
the sign of charge). Mass loading (inertia) or dissipation will reduce electric field to E < B. 
It is commonly assumed (mostly for the purposes of numerical simulations [9, 10, 22]) that 
the regions where E > B will be strongly dissipative due to resistivity. 

We can estimate the volume where E > B as « RqPo- Dissipation of the electromagnetic 
fields inside the regions E > B will lead to energy flux into those regions L diss R 2 g EqBq 
(this estimate neglects the fact that not all the surface of the black hole is covered by the 
regions with E > B, see Fig. 3). Thus, qualitatively, there will be local dissipation of the 
magnetic energy flux through the Schwarzschild circle. For merging black holes the resulting 
power 

L di ss*RlhBl = Bl^^ (44) 

For R 3> Rg this power is somewhat larger than the Poynting power of both the unipolar 
inductor and the Faraday wheel mechanisms (1), yet the estimates that went into Eq. (44) 
are likely to be solid upper limits due to the neglect of the geometry of the E > B regions, 
which, in fact, does not cover most of the black hole surface for ft <C 1. In addition, a 
fraction of this power will be swallowed by the black hole, while the escaping part will be 
heavily redshifted and unlikely to be observed. 



VI. FORCE-FREE MAGNETOSPHERES OF SCHWARZSCHILD BLACK HOLES 

We have argued above that the parallel electric fields created by the black hole motion 
lead to vacuum break down, generation of dense plasma that in turn screens the parallel 
component of the electric field. If the matter energy-density is much smaller that the energy- 
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density of the magnetic field, the plasma behavior will be controlled by magnetic field, while 
nearly massless charge carriers provide the currents demanded by the dynamics and ensure 
E B = condition. This is called the force-free limit [32]. Using the 3 + 1 formulation of the 
General Relativity [16], Eq. (7) , taking the total time derivative of the constraint E ■ B = 
and eliminating D t ~E and D t ~B using Maxwell equations, one arrives at the corresponding 
Ohm's law in Kerr metric 

. (B-Vx(aB)-E-Vx(«E))B + «(V-E)ExB 

J = 4^ (45) 

Note that this expression does not contain function /3. Eq. (45) is a non-linear equation 
for the time-dependent structure of magnetospheres of Kerr black holes. In the stationary, 
V x (qE) = 0, axisymmetric case Eq. (45) and Maxwell's equations (7) reduce to the 
Grad-Shafranov equation in Schwarzschild metric [15] for the A$ component of the vector 
potential. In our case, the lack of a cyclic variable prevents the reduction of Eq. (45) to a 
single equation. 



VII. DISCUSSION 



The electro-dynamics of Schwarzschild black holes moving through constant magnetic 
field resembles in many respect the pulsar magnetospheres. The motion of a black hole in 
vacuum both generates non-zero electric field along the magnetic field, and, in addition, 
produces regions where E > B. Similarly to the case of rotationally powered pulsars, the 
non-zero parallel electric field will lead to vacuum breakdown, generation of pair plasma and 
production of large-scale electric currents. Magnetospheres of moving black holes will have 
pair formation fronts and outer gaps, where the sign of the induced charge density changes. 
There is an important difference between pulsars and black hole: in case of pulsars the 
parallel electric field is produced by real surface charges [18], while in case of black holes the 
parallel electric field is a pure vacuum effect, resulting from the curvature of the space-time. 

In case of merging black holes we expect two kinds of electromagnetic outflows: those 
originating in the close vicinity of the black holes due to unipolar induction mechanism and 
those originating within the orbit due to the rotation of space-time. Both types of jets have 
approximately the same power. There is a clear thought experiment where the two effects 
(due to rotation of space-time and due to linear motion of a black hole) differ. Consider 
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a linear motion of a black hole through magnetic field. The effect discussed in this paper 
will appear, while the one discussed in Ref. [6] would not. On the other hand, consider 
a rotating massive ring: in this case the effect considered in Ref. [6] will appear, but one 
considered in this paper would not. 

The power of a unipolar inductor is taken from the energy of the linear motion of the BH. 
Thus, there is an effective friction force exerted by the magnetic field onto the black hole. 
Qualitatively, there are two ways to create a Poynting flux: (i) in vacuum due to changing 
electromagnetic fields; (ii) in plasma by generating a current-carrying plasma outlfow. In 
vacuum there is no energy loss by a black hole: even though the magnetic fields are disturbed 
by the passage of a black hole, there is no wave emission, as can be seen from the fact that 
in the frame of the black hole the electromagnetic fields are stationary. 

In the case of plasma, no time-dependence is necessary in order to produce a Poynting 
flux: presence of plasma allows one to chose from all the possible reference frames a one 
special frame, where the electric field is zero. Only in that special frame there is no Poynting 
flux, all other frames will have Poynting flux. Mathematically, in case of plasma the stress- 
energy tensor is diagonalizable, while this is generally not true in vacuum. Overall, the 
system under consideration is very similar to pulsars, in vacuum aligned rotator does not 
spin-down, but in reality parallel electric fields generate currents that carries electromagnetic 
energy, extracted from rotation. 

Any attempt to simulate numerically the magnetospheres of moving black holes using 
MHD-type (fluid) codes will face a problem of non-zero E-B ^ and the regions with E > B. 
Both these conditions violate a commonly used fluid assumptions. Pair production resulting 
from E-B^O will eventually ensure that E • B m in the bulk of the plasma. Thus, one 
possible way to simulate the black hole magnetospheres (again following the work on pulsar 
magnetospheres) is to assume that the regions where the ideal approximation is violated 
are sufficiently small, so that in the bulk the condition E • B = is satisfied. In the highly 
magnetized limit the dynamics then reduces to the force-free limit §VI. But experience in 
modeling pulsar magnetospheres [22] tells us that in the highly magnetized limit the system 
would evolve towards violations of the ideal condition E ■ B = by creating current sheets, 
where the magnetic field reverses, and, in addition, would spontaneously create regions with 
E > B. Resolving current sheets, or finding a proper prescription for treating them was 
a major obstacle in numerical simulations of pulsar magnetospheres. The appearance of 
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regions with E > B required introduction of artificial resistivity and demonstrating that the 
final outcome is largely independent of these ad hoc numerical procedures). This was the 
approach taken in Ref. [9-11] who performed a number of force-free simulations of black 
hole magnetospheres. Our results are generally in agreement with Refs [9-11]. 

Finally, we note that in the standard model of particle physics the non-zero second 
Poincare electromagnetic invariant leads to the appearance of sources of topological axial 
vector currents that can lead to the local violation of the baryon and lepton numbers through 
the triangle anomaly (which is responsible, e.g., for the two photon decay of it ). The 
triangle anomaly violates baryon number through a nonperturbative effect [33, 34]. In case 
of a black hole moving through a constant magnetic field, there is a non-zero divergence of 
the electromagnetic topological current J v 

Jv = A»(*F, U ) 
J = A • B = 

Ji = E x A + — B 

a 

7 M 7 

= -- sin26 cos <j)B E — = -E B (46) 

Note that the helicity, the time component of the topological current, is zero; the anomaly 
appears due to the 3-divergence of the spacial components of J v . We leave a more detailed 
investigation of the implications for the standard model of particle physics to a future work. 

I would like to thank Sergei Khlebnikov, Sergei Komissarov, Luis Lehner, Jonathan McK- 
inney for discussions. 
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Appendix A: Black hole magnetospheres in Kerr-Schild coordinates 

Interpretation of results in General relativity is often a non-trivial exercise. The results 
of the main part of the article were obtained in Schwarzschild coordinates. Schwarzschild 
coordinates have a singularity on the horizon, which often makes interpretation of the results 
problematic. To avoid the singularity the Kerr-Schild coordinates are often used instead. In 
this appendix we repeat the previous calculations in the Kerr-Schild coordinates and show 
that qualitatively the above-derived results generall hold. 

The Kerr-Schild metric is given by 

fa* = _ a 2 dt 2 + ^ dtdr + ^ + ^ dr 2 + r 2 ( d 02 + gin 2 ^ (M) 

The vacuum Maxwell equations (3) then gives equation for the same as in Schwarzschild 
coordinates, Eq. (4). Thus, in Kerr-Schild coordinates the flux function corresponding to 
constant field at infinity is A^= (B /2)r 2 sin 2 9. 

Using the covariant V operator with e r = , 1 d r we can then find the electromagnetic 
fields and the charge density 

E = \ j 1 + — VA = E [ cos 9'e r 1 sin 9'e e 

B = V x A^e^ = B \ cos 9e r . sin 9e e 

yl + 2M/r 

E M sing cos ((M + |) S jn 2 (g) + f + l) 



Qualitatively, the results in the Kerr-Schild metric look very similar to the ones in Schwarzschild 
metric. The only noticeable difference is the location of the E = B surface, which in the 
case of the Kerr-Schild metric is given by 

2M (sin 2 9 - ft (cos 2 9 cos 2 + sin 2 0) ) 



i - k (A3) 



see Fig. (8). 
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